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I. INTRODUCTION 

Quantum entanglement has emerged as a key re- 
source for quantum computing, quantum communication 
and information related processes. Characterization and 
quantification of multipartite entanglement is an inter- 
esting challenging problem in the field of quantum infor- 
mation and computation. There are several approaches 
to deal with the issue and various measures have been 
proposed (see e.g. [1], [2] for nice reviews). For mul- 
tipartite pure states, a frequently used measure is the 
Geometric measure of entanglement (GME) [3-8]. GME 
is a distance measure and for an arbitrary pure state |^), 
it is usually defined as 

GME{m = l-K% (1) 

where 

A^= max \{4>W)\ ( 2 ) 

\4>) e { product states } 

The notion of GME has also been extended to mixed 
states through the convex-roof construction. Another 
well-known measure, the Groverian measure of entangle- 
ment [9-11] is exactly the same as GME, up to a square 
operation. The Groverian measure originated from the 
modified quantum search algorithm [9] and is defined 
through the success probability in the search using the 
given state as the initial state: 

GM) = y/1 -PmaAm (3) 

where 

Pmax{\il>))= max |(</#)| 2 (4) 

\4>) e { product states } 

Through out this paper we shall consider the definition 
(3) for GME. 

It is clear that analytic calculation of GME for an ar- 
bitrary state is reasonably difficult since it involves non- 
linear maximization process. However, recently it has 
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been shown [5] that if a pure state is permutationally 
invariant, then the calculation of GME becomes greatly 
simplified. In this case, the optimal product state can 
be taken as tensor product of the same single system 
thereby drastically reducing the number of variables in 
the maximization process. 

Like all other well-known entanglement measures (von 
Neumann Entropy for bipartite pure states [12], Nega- 
tivity [13], Concurrence [14] etc.) there are some known 
bounds on the GME of superposition of pure states (see 
e.g., [15, 16]). In this work, we will explicitly calculate the 
analytic expression of GME for arbitrary superposition 
of two ./V-qubit canonical orthonormal GHZ states and 
also for W states. As we shall see later, these states do 
not saturate the bound obtained in [16]. Clearly, the for- 
mer states are not permutationally invariant in general 
and hence the theorem given in [5] can not be applied 
directly. However we will use it to derive our result. 

There exist quantum correlations other than entangle- 
ment, and one such is the quantum discord which has 
been studied extensively [17]. Basically, it quantifies the 
total non-classical correlations in a quantum state. How- 
ever, most of these investigations pertain to bipartite sys- 
tems only. Here, our aim is to study multipartite states 
and hence we shall follow the approach of Modi et. al. 
[18] in quantifying discord. In this unified view of all 
kinds of correlations in a quantum state, the quantifica- 
tions are done by the relative entropy, which makes this 
approach more challenging. Indeed optimization of rela- 
tive entropy is known to be a difficult problem, e.g., there 
does not exist a formula for the relative entropy of en- 
tanglement (REE) even for the simplest case of arbitrary 
2-qubit (mixed) states [19]. Very recently, the reverse 
problem of finding the set of entangled states for which 
the chosen separable state is the closest one, has been 
solved [20]. The necessary and sufficient condition given 
therein would help us in our analysis. 

The organization of this article is as follows: In Sec 

II, we will calculate the explicit expressions for GME of 
superposition of two orthonormal iV-qubit GHZ states 
and also of two W states and compare the results. In Sec 

III, quantum discord is derived analytically for the former 
class of states. In Sec IV, we conjecture the discord for 
iV-qubit W states to be log 2 N. Finally we discuss some 
possible extensions and related issues in Sec V. 
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II. GME OF ARBITRARY SUPERPOSITION OF 
TWO GHZ AND W STATES 

In this section we consider the superposition of two 
./V-qubit GHZ states and derive its GME. An analogous 
analysis is done for W states and the results are com- 
pared. 



A. GME of arbitrary superposition of two iV-qubit 
GHZ states 



Noting that the right hand side of (9) can be written as 

n ) + c^\<f) n }) , it follows 

that [22] 



max{ 



cos q sin 



2 ' 2 

^ if m = 1 or ?i = 1 



^} if m,n > 2 



(10) 



The only known bound for GME of superposition (Eq. 
(3) in [16]) gives P ma x < \ + cs, where c = cos a, s = 

sin a. Since for cs ^ 0, max{|-, 4-} < i, therefore P mox 
can not saturate this bound. 



The full set of canonical orthonormal iV-qubit GHZ 
states is given by (up to an irrelevant global phase) 



\B N (i))±\B N (2 N 
V2 



(5) 



where i = 0, 1, ... , 2 JV_1 — 1 and |.B/v(«)) = |h«2 • ■ -ijv) 
is the 'binary representation of the decimal number i in 
an iV-bit string'. In a more convenient notation (5) can 
be written as 



|0z 2 i 3 ■ ■ An) ± \U2i3 ■ ■ - »jv) 
V2 



(6) 



where ifc = 0, 1 Vfc = 2, 3, . . . , N and a bar over a bit 
value indicates its logical negation. 

Let us now consider an arbitrary superposition of two 
orthonormal GHZ states \Gf) and |G^) (i 7^ j) as follows 



|G y )=cosa|Gf) 



sin ae n \Gf) 



(7) 



First of all we note that the state in (7) is not Schmidt 
decomposable (SD) [21]. From the definition, if \ip) is 
SD, then 



y/Xi\ii 



Schmidt {A^}. 



(8) 



So, tracing out any subsystem would give a separable 
state and it is a necessary condition for SD. However, 
in this case, tracing out the third qubit from the 3-qubit 
state |Goi) = cosa|Gj)+sinae I7 |G^") gives an entangled 
state, which can be easily checked by the PPT criterion. 

We also note that (7) is not permutationally invariant 
and hence the theorem in [5] is not directly applicable to 
it. Though apparently it looks that the phase 7 can not 
be driven out, we will show that the GME is independent 
of 7. 

To calculate the GME, let us write the state (7) in the 
following convenient form 



ijl — Cl\4> m )\(j)n)+C2\4l m )\4>n)+C3\4> m )\<j)n)+C4\'llj m )\4>n) 



(9) 

where \tp m ) is the collection of (qu)bits when the two 
strings \Bn{i)) and |Bjv(j)) agree and \4> n ) is the collec- 
tion when they disagree, m + n = N; C\ = = ±c 4 , 

± e ' Ts i na = ±c 3 . 



B. Proof of the result (10) 

By suitable local unitary (LU) operations the state (9) 
can be transformed to the following form 



\Gi 



C-2 



V2 



-^(|0) TO |0)„±|l) m |l) n )+-^e < T(|0) m |l) n +|l) m |0) n ). 

(11) 

Clearly (11) remains invariant under permutation 
among the first m parties (and similarly for the rest n 
parties) and hence we can assume the optimal product 
state to be of the form [23] 

|$) = (cos#i|0)+e lAl sme 1 \l))® m (g>(cos6 2 \Q}+e lX2 sin0 2 |l»® n 

(12) 

with < 0i, 2 < f and < Ai, A 2 < it. 

In order to get P m ax{\Gij)), we have to maximize the 
quantity 

2Q = I cos™ 0i(ccos" 2 + se^-" A2 > sin™ 2 ) + 
e -imAi sin m 6 l 1 (±ce- mA2 sin™ 2 + se <7 cos™ 2 )| 2 (13) 

Now, applying \z± + z%\ < \z\\ + l^l successively, we 
have 

I cos™ 0i (c cos" 2 + se 4(7 -" A2) sin" 2 ) + 
e - mAl sin" 1 0i(±ce- mA2 sin" 2 + se^ cos" 2 )| 

< cos™ 1 |(ccos" 2 + seJ^- nX ^ sin™ 2 )| 
+ sin" 1 1 |(±ce- mA2 sin" 2 + se 47 cos" 2 )| 

< cos™ 0i (c cos" 2 + s sin" 2 ) 

+ sin'" 0i (c sin" 2 + s cos" 2 ) (14) 

Note that 

max [cos 2 0i (c cos 02 + s sin 02) 

O<01,02<f 

+ sin 2 0i(csin0 2 + scos0 2 )] = 1 (15) 

which occurs at (0i, 2 ) = (0, a), (f , f — a); additionally 
at ( arbitrary 6\, ?) if a = j. 

So taking n = 1 (the case m = 1 is similar by symme- 
try) in (14), 

cos" 1 0i (c cos 02 + s sin 2 ) + sin™ 0i (c sin 2 + s cos 2 ) 

< cos 2 0i (c cos 02 + s sin 02 ) + sin 2 0i (c sin 02 + s cos 02 ) 

< 1 (16) 
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where, in the second line, we have used cos™ 9 < 
cos 2 9, sin m 9 < sin 2 9 Vto > 2. Thus, if n = 1 or 
m = 1, we have from (13), (15) and (16), P m ax(\Gij)) = 



Similarly using 

max [cos 2 6i (c cos 2 92 + s sin 2 #2) + sin 2 9\ 

O<01,02<f L 

(csin 2 02 + s cos 2 #2)] =max{c, s}, (17) 
it follows that P TO ax(|Gy)) = maxfy, 4^-}, if m, n > 2M 



C. GME of arbitrary superposition of two iV-qubit 
W states 

The GME of superposition of two 3-qubit W states has 
been presented in [3] . Although the generalization to N- 
qubit case is quite straightforward, still we wish to derive 
it for completeness and it will help us to compare with 
the case of GHZ states. [Interestingly, it is shown that 
the 4-qubit case is the easiest one (easier than even the 
3-qubit case) and an explicit formula of GME has been 
derived] . 

Let us consider an arbitrary superposition of the two 
iV-qubit W states namely \W) = ^=(|00 . . . 1) + . . . + 

|10...0)) and \W) = ^j-(|01...1) + . . . + |11 . . . 0)) as 
follows 



\WW) = cosa|W0 + e 17 sina|W) 



(18) 



We note that the LU transformation {|0),|1)} — ^ 

{|0), e~ N - 2 11)} leads to an overall phase and hence the 
phase 7 becomes redundant. So without loss of general- 
ity, we can assume 7 = and therefore 

\WW) = cos a\W) + sin a\W), 0<a<^ (19) 

Since the state |WM^) is permutationally invariant 
(with positive coefficients), we can assume the nearest 
product state as [5] 

|0) = (cos 6\0) +sin0|l»® J>r J O<0<! (20) 

The overlap is given by 

WW)\ 2 = Q = NlsC^S + cCS N - 1 } 2 (21) 

where C = cos 9, S = sin (9, c = cos a, s = sin a. 
The condition for maximum of Q (^ = 0) becomes 

s[t N - (N - l)t N - 2 ] + c[(N - l)t 2 - 1] = (22) 
and hence we have 



PmaxQWW)) 



Nt 2 



(l + t 2 ) N 



[c + st 1 "- 2 ] 2 (23) 



where t = tan# (> 0) will be determined from the poly- 
nomial equation (22). We note that Eq.(22) has only one 
positive root for TV = 3, 4 and has at most three positive 
roots for N > 5. Hence the GME can be calculated using 
numerical techniques of root finding. 
Particularly, for N = 4, we have 



t 2 



3(c 



V9 - 14cs 



2c 



(24) 



which readily gives the expression for GME. The graph 
of GME vs. s for 4-qubit case is exactly analogous to 
the 3-qubit case ( Fig.-l. in [3]) and hence we are not 
reproducing it here. 

It is worth mentioning that if we consider the super- 
position of two other (even LU equivalent) W states, the 
GME will not be the same. 



D. Comparison between GME of superposition of 
GHZ and W states (with respect to the results 
presented here) 

1. For N = 3, either m = 1, or n = 1 and hence 



Pmax ( I ) ) 



,(\Gt)) 



(25) 



Thus for the three-qubit GHZ states, the GME of 
the superposition is independent of the superposi- 
tion parameter a and the phase 7 whereas for the 
W states, it is dependent on a. Of course, the 
GME of superposition of \Gf) and |G~) (which is 

mill{ | co S a+mnae^ |J cos a -sin a£Z |} ) depends Oil both a 

and 7. 

For arbitrary TV, if m = 1 or n = 1 (i.e., if the 
Hamming distance between \Bn(i)} and \Bpf(J)) be 1 
or IV - 1), then G(|G«)) = ^. 

2. For N = 3, we note that by superposing two 
orthonormal W states, we can get the resultant 
entanglement equal to that of a GHZ state. For 
■ \ w ^)+\ w -) \ J- G(\Gf}), where 



example, G (™^) = ^ 



\Wi) 



±=(|001) - w|010) + w 2 |100)),|W 2 



^(-|001) + w 2 |010) - w|100)); to being a complex 
cubic root of unity. If we consider the superposition of 
W and W, it follows from Fig-1 in [3] that we can choose 
a specific value of s to get G{\WW)) = ^ = G{\GHZ)). 
On the other hand we can not get the entanglement 
of a W state by superposing any two orthonor- 
mal GHZ states from the canonical set (5) since 
G(|G«» < ^ < G(\W)). 

However, for N = 4 the situation is different. 
By superposing two W states, we can get the en- 
tanglement of a GHZ state and vice-versa. For 

GQGf}), where 



example, G( 1MM ) 
\Wi) 



1 

72 



1(10001) + |0010) + 10100) + 1 1000)) and 
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\w 2 ) = Kloooi) + |0010) - |0100) - 1 1000)). It 
can be checked that by superposing |G ) and 
IG3) we can get the entanglement of the W state: 

& = G(\W)), where |G 03 ) 



with a given by sin a (or cos a) 



is given by (7) 



For large N, the situation is somehow opposite to the 
3-qubit case. Here we can always get the entanglement of 
a W state by superposing two GHZ state, but we don't 
know if the converse is also true. This requires further 
investigation. 

3. For JV = 3, even if we consider equal superposi- 
tion (a = ?)7 = 0), the state \Gij) is not invariant un- 
der permutation. The corresponding product state would 
also not be permutationally invariant in general. How- 
ever, for some \Gij), the optimal product state may still 
be permutationally invariant. As an example, the state 
IG&) = |(|000) +]lll) + |001) + [110)), is not permuta- 
tionally invariant. But we can choose the nearest product 



state as 



- |Q)+|i) \®3 



V2 



)® 3 . [Of course, there exist other optimal 
product states which are not permutationally invariant 
e.g. \q)\q)\+) where q = 0,1,- and |±) = 
This may lead to the intuition that the optimal prod- 
uct state may be permutationally invariant even if the 
state itself is not. However, in [5], the authors have 
proved the stronger result that this is not the case if 
the state is genuinely entangled (or we look into those 
parties for which this state is entangled and symmetric, 
e.g., |Gqi) = |</> + )|+), so we have to consider only the 
first two particles. Being symmetric and entangled, the 
combined state of these two particles have the necessarily 
symmetric closest product states |s)|<?))- 



III. QUANTUM DISCORD FOR THE CLASS OF 
STATES (7) 

As mentioned earlier, we will follow the approach of 
[18] to characterize and quantify all kinds of correlations 
in a quantum state. The definitions of relevant quantities 
are: 

Entanglement E — min S(p\\a) (26) 

trE'D 

Discord D = mm S(p\\x) (27) 

Dissonance Q = mm S(a\\x) (28) 

Classical correlations C = min S(p\\tt) (29) 

ttEV 

where V is the set of all product states (i.e., states of the 
form 7r = 7Ti <£> 7T2 <g) . . . <E) itn ) , C is the set of all classical 
states (i.e., states of the form x = J2~£ P~g\ k)(k |, with 
the local states \k n ) spanning an orthonormal basis), T> 
is the set of all separable states (i.e., states of the form 
a = J2kPk n i ® ■ • -®^n) and s ( x \\y) = Tr[iTogx- 
x log y] is the relative entropy of x with respect to y. We 



shall first find out the closest separable state (CSS) to 
the class of states (7). Fortunately it turns out that the 
CSS is also a classical state thereby implying D = E = 
ReQGh)), Q = 0. 

Before proceeding to calculations, we recall that find- 
ing out the CSS is a challenging problem [19, 24]. To 
obtain the CSS to a multipartite state, two interesting 
tools are available in the literature. The first one is a 
lower bound through the generalization of Plcnio-Vcdral 
formula [25]: 

S(p N \\a N ) > S(pjv_i||(Tjv_i) + S{p N -i) - S(p N ), (30) 

where pm is any Af-partite state and er/v is an ^-separable 
state. So, for any Af-qubit pure state pn we have the 
lower bound 



E(p N ) > max{ J B(p Ar _ 1 ) + 5(pat_i)} 



(31) 



where the maximum is taken over all possible bipartition 
of N — 1 versus single qubit. 

The second tool is due to Wei et. al. [4]. For any 
A^-partite pure state it gives a lower bound on E 
through P max (\tp)): 



E(\tP)(iP\) >-log 2 P mox (h/>)). 



(32) 



Since E is defined through minimization, if we can find 
a separable state a which saturates either of the bounds 
in (31) and (32), then a will be the required CSS. In- 
fact the later bound has been extensively used to de- 
rive REE of symmetric Dicke states [4] and even mix- 
tures of them [26]. But unfortunately, these bounds are 
not saturated for the states in (7). Indeed, the bound 
(32) is not saturated even for the simplest 2-qubit non- 
maximally Bell states (e.g., for \<j>) = a|00) + 6 1 1 1) with 



E{\4>)) 



\b\ 2 = 1, one has P„ 



max{|a| 2 , |6| 2 } whereas 



; log|a| 2 - H 2 log|6| 2 = H{\a\ 2 ) [27]). It 



is thus quite challenging to derive the CSS. The reverse 
problem (i.e., starting from a a on the boundary of T>, 
determining all entangled state p for which a is the CSS) 
is also interesting and has been solved for the 2-qubit 
case [28], very recently for multiparty states [20]. We 
shall apply this multi-party criteria to derive the CSS. 
The criteria reads: 

Necessary and sufficient criteria for CSS [20]: 
a 6 T> is a CSS for an entangled state p if and only if 



max Tr a' ' L a {p) = 1, 

u' ED 



(33) 



where the linear operator L a is defined in the following 
way. Let the eigendecomposition of hermitian positive 
operator a be a = diag(ai, aa, . . . , a n ). Then for any 
P=K]i, j= i, L a (/3) is defined by 



bklT;, 



if dfc = en = a 



(34) 



We shall now derive the CSS of our states. Since REE 
is invariant under LU and (7) can be transformed to (11) 
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by LU, we can consider REE of this state, without loss 
of generality. The state (11) has GME similar to the 
non-maximal Bell state. So we assume that it will have 
a similar REE also. Hence we take the CSS as 

o = |(|00)(00|+|11)(11|)+^(|01)(01| + |10)(10|), (35) 

where (and hereafter) we have dropped the suffixes m 
and n. 
Proof: 



a = diag( 



c 2 ,s 2 s 2 c 2 



2 ' 2 

^,2 



and p = — 
' 2 



2 ' 
cse" 



cse 



-»7 



cse 



vy 



\ 



±C 2 \ 

±cse 11 

^,2 



/ 1 






±1 




qe l ~< 


1 


1 


±qeP 




qe l ~< 


1 


1 


±qe i ~< 




V ±1 


±qe- l ~< 




1 


) 



±c 2 ±cse~^ ±cse~ i 'i 
Hence from the definition of L CT (/5), 



lap) 



cs log —tj . . 

where q = — i _ § . Note that \q\ < 1. 
Now let a' = Y,Pk\<t>k){4>k\- Then 

Tr o'L a {p) = ^2p k ((j> k \L a (p)\(j) k ) 

= ^ Pfc [|(0 fe |OO)| 2 + |(0 fc |Ol)| 2 + |(0 fe |lO)| 2 + |(0 fe |ll)| 2 

+ 2Real{qe~^{4>kmm\4>k) + (lO|0 fe >)) 
± 2i2cal<^ fc |OO><ll|0 fc > + 2i2eoJ<0 fc |Ol><lO|0 fc > 
± 2fea/(g e -^(^|ll)((01|^) + (lO|0 fe )))] 

< ^p fc [|(0 fe |OO)| 2 + |(0 fc |Ol)| 2 + |(0 fe |lO)| 2 + |(<fe|ll)| 2 

+ (|(^|00)||(^|01)| + ... + |(^|10)||(^|ll)|)] 

= 5> fc [|<&|00>| + |(^|01)| + |(^|10)| + |(0 fe |ll)|] 2 

Since each \4> k ) is a product state, we have \4> k ) = 
l^fc) |"0fc) - So the last expression above can be written as 

E«b[(Kwb|o>| + K^|i)l)(K^|o)| + \mi)\)f < i, 

since for any normalized product state \<j>) (of > 2- 
qubits), |(0|O)| + |(0|1)| < 1 (which can be seen from 
(17)). ■ 
Thus a is indeed the CSS. Being a classical state as 
well, a is also the CCS, thereby yielding D = E = 

-c 2 log ^ - (1 - c 2 ) log ( -^p- = 1 + H{c 2 ) and Q = 0. 
We have depicted all the known bounds and our exact 
results for this state in Fig-1. 



IV. CONJECTURE FOR DISCORD OF 
iV-QUBIT W STATE 

From the discussion of the previous section, it is clear 
that determining CSS is a non trivial task. Determining 




FIG. 1. (color online only) Known bounds Vs. exact results: 
BG = bound on GME obtained in [16], G = exact GME, BE= 
bound on REE obtained through — log 2 P ma x [4] , E = exact 
R EE 



the CCS is even more complicated because the set C is 
not a convex set and hence the standard tools of convex 
optimization theory is not directly applicable. However, 
to calculate the discord D and dissonance Q, the au- 
thors of [18] have simplified the task of minimizing over 
C. They have shown that for any given p, the CCS \ p is 

given by x P = I k K k \p\ k X k l> where {| k )} forms 
the eigenbasis of \p- This simplifies expressions for D 
and Q as the minimization of the relative entropy over 
C reduces to minimization of the von Neumann entropy 
S(xx) over the choice of local basis {| k )}: 

D = S(x P )-S(p), Q = S( X *) - S(a), (36) 



where S(xx) = min 



S[\t)(t\x\t}(t\). Therefore, 



for numerical computation of D, one can choose arbitrary 
local bases and find the minimum of the corresponding 
entropies. An even finer approach is to generate a vector 
(with equal spacing) and using Gram-Schmidt method 
construct a complete orthonormal basis and obtain the 
minimum entropy. This technique is useful mostly in low 
dimensional cases [29]. 

The CSS to the A-qubit VK-statc is known to be [4] 



JV 



fe=0 



fc\ fe fN -k s 



A 



A 



\S(N,k))(S(N,k)\, 



\S(N,k)) being the fc-th symmetric (Dicke) state. For 
A > 3, the above separable state is not a classical state. 
Therefore D ^ E and Q ^ for W states (contrary to 
the GHZ case, where the CSS was a classical state). 

Since the W state is symmetric, we assume that the 
CCS can be chosen to be symmetric [30]. So we choose 
each of the local orthonormal basis of the classical state 
Xw as 



|C) 
II') 



Vp|o> + V^pW 
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FIG. 2. (color online only) p Vs. 



so that (x'\x) = {-l) x ^p, (x'\y) = s/T=p, x^y = Q,l. 
Therefore we have 



. ..x' N \ Viy2 ... VN ) = (-lr^nyw) 



N- 



where m is the number of positions where the two binary 
strings x and y agree and m\ is the number of positions 
where they both have 1. Since for the W state each y has 
exactly one 1, the inner product (x'|VF) will just depend 
on the number of l's in x. So, if a basis \xk) has k number 
of Is, we have 



i 

N 

+ 



N - k c 1 ( y /p) lf - k -Hy/T=p) k+1 



= —(^-*-\jr^,f-\N{l -p) k] (37) 



Now from (36), to determine D, we have to find the 
minimum of 



53 \x')(x'\W)(W\x')(x'\ 
J2 \(x'\W)\ 2 \og 2 \{x'\W)\ 2 

x 2 ---% N — 0,1 

^ N C k X k \og. 2 X k , (38) 



TV 



fc=0 



where = Kx'jVF)! 2 with x k being any binary string 
of length N having k Is, is given by (37). It can easily 
be checked that S has (global) minimum at p = 0, 1 (see 
Fig.-2). Therefore the CCS to W state is the dephased 



state in computational basis and consequently we have 
D = log 2 N. 

Employing the method of [29], we have also numer- 
ically verified (independent of the assumption that the 
CCS is symmetric) that upto N = 5, this indeed is 
the minimum. We thus conjecture that discord of 
N-qubit W state is log 2 N. 

V. DISCUSSION 

First of all, we note that both the results (10) 
and (35) can straightforwardly be extended to the 
case of non-maximal GHZ states (i.e., a\iii2 ■ ■ ■ ijv) + 
b\iii2 ■ ■ ■ In), |a| 2 + \b\ 2 = 1). However, calculation of 
GME for superposition of two arbitrary GHZ states is 
more involved. In fact, even for a single non-maximal 
(generalized) W state, obtaining the GME is quite non- 
trivial. Recently, the three qubit case has been studied 
in [6] which has been further generalized to iV-qubits [7] . 
From a broader perspective, a generalization of GME in 
which the maximum distance would be calculated from 
the set of all states which arc equivalent under stochastic 
local operations and classical communications (instead 
of just product states), has been introduced in [31]. It 
would be interesting to see how the GME of superposi- 
tion behaves in this context. 

Another basic question related to the measure of cor- 
relations is the additivity of the proposed measure. It is 
known that GME is in general not additive [32] ; precisely, 
for N > 3, GME is not additive for any two iV-partite an- 
tisymmetric states [33]. However, this is still not known 
for the case of total correlations T p ( defined as S(p\\ir p )) 
in a quantum state. It has been conjectured [18] that T p 
is subadditive: T p > E+Q + C a , where C CT is the classical 
correlation S(x<t\\' k <t)- 

A further direction along our line of study would be to 
explore the correlations in TV-qubit GH Z-diagonal states 
(an arbitrary mixture of the states \Gf^){Gf-\). Because 
of the simple structure (both algebraic and geometric), 
the two qubit case allows easy computation of all the 
measures and has been studied extensively. But, beyond 
this, even the criteria for entanglement is unknown till 
date. We hope that the lower bound in (31) may pro- 
vide some insight in determining the structure of the CSS 
which then can be verified using the necessary and suffi- 
cient condition given in [20]. 

To conclude, we have derived analytically the GME 
and discord (via REE) for superposition of some or- 
thonormal GHZ states. We have also conjectured the 
discord for W states. Perhaps a similar approach could 
be applied to other permutationally invariant states. 
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